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1 .  INTRODUCTION 


The  purpose  of  this  effort  is  to  provide  quantitative  information  to 
aid  in  evaluating  the  effects  of  x-ray  spectral  variations  on  the  SGEMP 
response  of  satellite  systems.  The  problem  is  approached  by  defining  a 
number  of  excitation  parameters  such  as  electron  yield  and  velocity,  dose, 
dipole  moment,  etc.,  which  are  the  driving  functions  for  the  various  SGEMP 
effects.  The  driving  functions  are  evaluated  at  various  locations  in  a 
representative  satellite  structure  ranging  from  the  external  surface  through 
various  structural  elements  to  heavily  shielded  equipment  boxes. 

The  excitation  parameters  chosen  are  basically  linear  and  can,  there¬ 
fore,  be  computed  as  a  function  of  incident  x-ray  energy.  Coefficients  are 
computed  which  relate  the  various  excitation  parameters  to  incident  x  rays 
such  that  the  excitation  parameters  can  be  obtained  for  an  arbitrary  spec¬ 
trum.  The  effects  of  various  x-ray  spectra  can  be  evaluated  by  comparing 
the  excitation  parameters  produced  by  candidate  sources  with  those  produced 
by  blackbodies  or  weapon  spectra  of  interest.  All  excitation  parameters  are 
computed  for  blackbody  spectra  ranging  from  1  to  15  keV. 

While  most  of  the  excitation  parameters  chosen  are  linear  with  x-ray 
fluence,  it  should  be  noted  that  many  of  the  x-ray  responses  are  nonlinear. 

It  is  of  some  interest  to  evaluate  certain  "response"  parameters  to  deter¬ 
mine  the  range  in  which  responses  become  nonlinear.  However,  the  nonlinear 
response  parameters  are  not  the  most  suitable  for  evaluating  x-ray  spectra 
because  the  effects  tend  to  become  less  sensitive  to  spectral  variations  in 
the  nonlinear  regime  (Ref.  1) .  Because  the  responses  are  more  sensitive  to 
spectral  variations  in  the  linear  regime,  and  because  linear  responses  occur 
at  numerous  locations  on  spacecraft  over  virtually  the  entire  range  of  inter¬ 
est,  the  linear  response  parameters  give  a  better  indication  of  the  maximum 
response  variations  caused  by  spectral  variations. 

The  nonlinear  responses  can  be  somewhat  artificially  divided  into  two 
categories  of  problems  —  the  internal  response  of  internal  cavities  and  the 
external  response  of  convex  objects  where  a  certain  fraction  of  the  electrons 
can  escape  from  the  object.  Spectral  effects  on  the  nonlinear  internal  cav¬ 
ity  response  are  evaluated  with  a  quasistatic  1-D  solution  with  emission 
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from  one  face  of  a  diode.  The  potential,  transmitted  current,  and  dipole 
moment  are  used  as  response  parameters.  The  external  response  evaluation 
in  the  nonlinear  regime  is  not  treated  here. 

Because  of  the  potentially  tedious  calculations  required  to  evaluate 
an  arbitrary  spectrum  and  compare  excitation  parameters  with  those  of 
blackbody  spectra  or  weapons,  a  set  of  magnetic  cards  for  the  TI-59  cal¬ 
culator  has  been  prepared.  These  cards  are  available  and  can  readily  be 
used  to  compute  each  of  the  parameters  considered  in  this  effort. 
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2.  SPECTRUM  EVALUATION  PARAMETERS 


There  exists  a  myriad  of  parameters  that  could  be  chosen  to  evaluate 
SGEMP  response  sensitivities  to  x-ray  spectrum  variations.  One  approach 
consists  of  identifying  all  important  known  effects  and  defining  responses 
(e.g.,  fields,  currents,  potentials,  etc.)  for  representative  configurations 
or  geometries  of  possible  interest.  Such  an  approach  tends  to  depend  upon 
the  configuration  or  model  chosen  and  might  also  result  in  an  excessive  list 
of  evaluation  parameters . 

A  second  approach  is  to  define  the  quantities  responsible  for  producing 
each  of  the  known  effects.  Presumably,  if  each  of  the  driving  functions  or 
excitation  parameters  is  specified,  the  responses  resulting  from  these  exci¬ 
tation  sources  will  also  be  uniquely  specified.  Implicit  in  this  treatment 
is  the  assumption  that  an  x-ray  spectrum  which  produces  excitation  terms  com¬ 
parable  to  a  second  spectrum  will  also  produce  similar  SGEMP  responses.  This 
approach  precludes  the  need  for  defining  complex  geometries  and  performing 
complicated  time-dependent  calculations,  and  was  chosen  for  this  effort. 

A  distinction  between  excitation  and  response  terms  may  be  helpful  at 
this  point.  Parameters  which  are  linear  with  fluence  and  independent  of  x-ray 
pulse  time  history  are  termed  "excitation"  parameters.  Parameters  that  become 
nonlinear  as  space-charge-limiting  becomes  appreciable,  and  are  dependent  on 
pulse  time  history,  are  called  "response"  parameters.  Furthermore,  source 
terms  are  dependent  only  upon  attenuation  and  emission  materials,  whereas 
response  terms  are  also  dependent  on  geometry . 

Primary  emphasis  in  this  effort  is  on  defining  a  pertinent  set  of  exci¬ 
tation  parameters  and  evaluating  them  at  suitable  points  on  a  representative 
spacecraft.  A  limited  set  of  response  parameters  is  considered  to  obtain  a 
first-order  characterization  of  spectral  response  to  space-charge-limiting. 
Solutions  are  kept  as  general  as  possible,  and  lend  themselves  to  analytic 
treatment. 

Variations  in  x-ray  time  history  are  not  specifically  considered  in 
this  exercise.  Effects  due  to  different  pulse  widths  of  spectra  can  be 
readily  estimated  in  the  linear  regime.  Evaluations  of  simultaneous  spectral 
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and  temporal  variations  in  the  nonlinear  regime  are  very  complex  and  beyond 
the  scope  of  this  simple  treatment. 

The  parameters  chosen  for  evaluating  spectra  are  summarized  in  Table  1. 
The  list  is  not  exhausive  but  is  representative  of  the  terms  thought  to  be 
important  in  producing  x-ray  responses.  The  parameters  can  be  conceptually 
divided  into  four  categories:  the  electron  yield  which  is  the  basic  x-ray 
interaction  with  the  materials;'  the  dipole  moment  in  a  dielectric  which  is 
descriptive  of  the  charge  deposition  profile;  the  energy  deposition  in  diel¬ 
ectrics  and  semiconductors;  and  finally,  several  quantities  such  as  average 
velocity  and  energy  and  dipole  moment  which  are  descriptive  of  the  location 
and  motion  of  charge  in  a  vacuum. 

Table  1.  Parameters  Employed  in  Comparing  SGEMP  Effects 
of  Different  Photon  Spectra 

Electron  Yield 
Primary  yield 
Forward-to-reverse  ratio 
High-to-low-Z  ratio 
Dielectric  Charge  Profile 
Dipole  moment 

Dose 

Bulk  dose 
Interface  dose 
Dose  ratio  (Au/Si) 

Vacuum  Charge  Profile 

Mean  electron  velocity 
Average  electron  energy 
Diode  potential 
Diode  current 
Dipole  moment 


In  terms  of  our  definitions  of  excitation  and  response  terms,  we  may 
categorize  electron  yield,  dose,  mean  velocity  and  energy,  and  dielectric 
dipole  moment  as  excitation  terms.  These  quantities  are  not  dependent  on 
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flux  rates,  and  are  linear  with  fluence  over  the  ranges  of  interest  here. 
Diode  potential,  current,  and  dipole  moment  are  classified  as  response 
terms.  These  require  a  specification  of  photon  source  flux  and  some  min¬ 
imal  system  geometry,  and  also  may  become  nonlinear  with  fluence. 

Four  basic  categories  of  x-ray  response  are  being  described  by  the 
parameters  of  Table  1.  They  are  the  external  SGEMP  response,  internal  cav¬ 
ity  SGEMP  response,  shielded  cable  excitation,  and  TREE  effects  caused  by 
energy  deposition  in  electronics,  as  indicated  in  Table  2.  The  electron 
yield  is  important  in  every  category  in  that  it  provides  the  basic  coupling 
of  photons  into  the  system.  Primary  electron  yield,  as  defined  here,  is 
the  total  number  of  electrons  produced  by  an  irradiated  surface  in  the 
energy  range  above  100  eV  (generally  speaking).  This  yield  is  distinguished 
from  secondary-electron  yield  (electrons  of  energy  100  eV  and  below)  because 
analytic  multiple-scattering  cross  sections  are  considered  to  be  inaccurate 
below  1  keV  (Ref.  2) .  Recent  analytical  work  in  the  regime  below  1  keV  (Ref. 
3)  has  not  been  incorporated  here. 

The  external  response  in  the  linear  regime  is  determined  by  the  primary- 
electron  yield  in  the  reverse  direction  and  by  the  average  electron  velocity. 
Under  certain  circumstances,  forward-directed  emission  may  be  important  also. 
Generally  speaking,  electron  emission  will  occur  from  low-Z  materials  on  the 
outside.  Stored  charge  release  of  abundant  low-energy  electrons  may  also  be 
an  important  SGEMP  driver  in  the  low-fluence  regime.  This  secondary-electron 
yield  has  been  shown  to  be  proportional  to  interface  dose  (Ref.  4) .  The 
spectrum  energy  becomes  important  in  the  nonlinear  regime,  but  the  response 
in  this  regime  is  less  sensitive  to  spectrum  than  in  the  linear  regime  and 
has  not  been  treated. 

The  internal  cavity  response  is  highly  dependent  on  electron  yield 
ratios  as  well  as  on  electron  velocity.  Cavities  may  be  large  (^2  m)  with 
emission  from  either  low-Z  or  high-Z  materials..  The  net  forward- to -backward 
emission  ratio  may  be  important  in  determining  net  current  in  cavities.  The 
larger  cavities  may  produce  considerable  space -charge -limiting  at  the  higher 
fluence  levels,  and  the  average  energy,  diode  potential  and  current,  and 
dipole  moment  may  be  used  to  characterize  the  response  (Ref.  5) .  Small  cav¬ 
ities  such  as  equipment  boxes  may  contain  higher-Z  materials  in  sizeable 
proportions.  High-to-low-Z  electron  emission  ratios  may  be  important  in 
these  cases. 
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Table  2.  Summary  of  Parameters  Describing 
Each  SGEMP  Regime  Treated  Here 


External  Response 
Linear  regime 

Electron  >ield 
Average  velocity 
Interface  dose 
Internal  Cavity  Response 

Large  cavities,  linear  regime 

Electron  yield  from  low-Z  materials 
Forward-to-reverse  ratio 
Average  velocity 
Interface  dose 

Large  cavities,  nonlinear  regime 
Diode  potential,  current 
Dipole  moment 
Average  energy 
Small  cavities 

High-to-low-Z  emission  ratio 
Cable  Excitation 
Electron  yield 
Dipole  moment 
Dose 

Energy  Deposition  in  Electronics 
Bulk  dose 
Interface  dose 
Dose  ratio  (Au/Si) 


The  excitation  of  shielded  cables  by  direct  x-ray  penetration  is  gen¬ 
erally  proportional  to  electron  yield  from  cable  shields  for  cables  with 
large  gaps,  and  proportional  to  the  electron  dipole  moment  resulting  from 
penetration  of  electrons  emitted  from  shielding  material  into  the  dielec¬ 
trics  for  cables  without  gaps  (Refs.  6,7,8).  At  the  higher  fluence  levels, 
dielectric  relaxation  times  may  become  comparable  to  response  times  of 
interest.  These  have  been  shown  to  depend  upon  dose  and  dose  rate  (Ref.  9). 

Transient  radiation  effects  in  electronics  (TREE)  are  dependent  upon 
energy  deposition  due  to  photo-Compton  electrons.  The  process  is  important 
both  in  the  bulk  (greater  than  an  electron  range  from  the  surface)  and  at 
material  interfaces.  In  the  latter  case,  dissimilar  materials  can  greatly 
affect  energy  deposition  within  an  electron  range  of  the  interface  (Ref. 

10) .  This  effect  is  represented  well  to  first  order  by  the  ratio  of  the 
bulk  dose  in  each  material . 

SGEMP  parameters  considered  in  this  exercise  but  not  directly  employed 
in  the  spectral  sensitivity  studies  are  listed  in  Table  3.  Certain  of  these 
quantities  are  not  strong  functions  of  photon  or  electron  energy  over  the 
range  of  1  to  1000  keV.  For  example,  electron  backscatter  yield  is  a  weak 
function  of  energy  spectrum  in  this  range  (Ref.  11).  Other  quantities  of 
interest  may  be  well  represented  by  a  different  parameter  type  as  far  as 
spectral  sensitivity  is  concerned.  Interface  dose  enhancement  due  to  dis¬ 
similar  materials  falls  into  this  category.  The  energy  sensitivity  of  the 
interface-to-bulk  dose  ratio  in  a  medium  behaves  similarly  to  the  bulk  dose 
ratio  of  the  two  mediums.  Therefore,  only  the  latter  ratio  has  been  included 
in  the  list.  Secondary-electron  yield  has  been  shown  to  be  proportional  to 
surface  dose  [energy  deposition  with  an  electron  range  of  the  surface  (Ref. 
4)],  and  so  has  not  been  computed  separately. 

Table  3.  Potentially  Interesting  Spectrum  Evaluation  Parameters 
Not  Treated  Here  Due  Either  to  Weak  Energy  Dependence  or 
to  Similar  Sensitivity  to  These  Quantities  Listed  in  Table  1 

Secondary-electron  emission  efficiency 
Electron  backscatter  yield 

Interface  dose  enhancement  due  to  dissimilar  materials 
Electron  angular  distribution 
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3.  SATELLITE  LOCATIONS 


Representative  satellite  materials  and  thicknesses  have  been  chosen 
with  assistance  from  TRW  (Ref.  12),  as  shown  in  Table  4.  The  intervening 
materials  are  listed  as  well  as  the  energy  at  which  the  photon  attenuation 
is  l/e.;  Also  indicated  is  the  blackbody  temperature  for  which  the  attenua¬ 
tion  is  1/e.  The  locations  have  been  chosen  to  range  from  the  very  thinnest 
(thermal  blanket)  to  the  very  thickest  (heavily  shielded  box),  with  a  good 
sample  of  intervening  thicknesses. 

The  materials  on  the  external  surface  of  the  satellite  tend  to  be  low-Z 
such  as  quarts  coverglass,  fiberglass,  and  aluminum.  Large  internal  cavity 
walls  are  also  characterized  by  low-2  materials,  although  a  significant 
amount  of  lead-tin  solder  is  found  on  solar  cells.  The  Faraday  cage,  thermal 
blanket,  and  solar  cell  models  represent  these  types  of  wall  materials.  The 
thin  mylar  or  thinnest  Faraday  cage  (3  mils  aluminum)  is  apt  to  produce  the 
largest  internal  cavity  response  due  to  maximum  electron  emission.  The 
thicker  Faraday  cage  model  and  solar  cell  give  a  measure  of  the  hardness  of 
a  spectrum  required  to  penetrate  thick  wall  materials. 

A  minimum- thickness  equipment  box  wall  (20  mils  aluminum)  is  also  mod¬ 
eled  to  produce  the  worst-case  response  driver  in  regions  where  vulnerable 
electronics  are  found.  Various  emission  and  dose  ratios  are  obtained  for 
this  minimal  attenuation.  Addition  of  10  mils  Kovar  represents  a  typical 
attenuation  through  the  semiconductor  case  into  the  interior,  and  maximum 
spectrum  attenuation  is  represented  by  the  heavily  shielded  box  model. 

Minimum  photon  attenuations  into  cable  interiors  are  represented  by  a 
shield  model  consisting  cf  3  mils  copper;  20  mils  copper  is  used  to  repre¬ 
sent  penetration  by  photons  to  a  cable  deep  in  the  interior  of  a  cable 
bundle. 
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Table  4.  Satellite  Locations 


Blackbody 

Energy 

Temperature 

[e-fold 

[e-fold 

Satellite 

Intervening 

attenuation] 

attenuation] 

Location 

Material 

(keV) 

(keV) 

Surface 

- 

Thermal 

blanket 

10  mils  mylar 

Faraday 

cage 

3  mils  aluminum 
(minimum) 

20  mils  aluminum 
(maximum) 

Solar  cell 

5  mils  quartz 

10  mils  silicon 

1  mil  lead-tin  solder 

5  mils  fiberglass 

2  mils  aluminum 

5  mils  fiberglass 

Cable 

shield 

3  mils  copper 
(min.  shield  thickness) 
20  mils  copper 
(behind  bundle) 

Box  wall 

20  mils  aluminum 
(minimum) 

Box  §  semi¬ 
conductor 
device  housing 

20  mils  aluminum 

10  mils  Kovar 

Heavily 
shielded  box 

20  mils  aluminum 

10  mils  Kovar 

12  mils  tantalum 
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4.  EXCITATION  PARAMETER  CALCULATIONS 


Satellite  excitation  parameters  have  been  evaluated  for  selected  sat¬ 
ellite  positions  and  materials.  Unit  fluence  is  assumed  to  be  incident  for 
all  configurations.  Results  are  valid  at  all  fluence  levels,  including  the 
nonlinear  response  regime.  A  modified  version  of  the  QUICKE2  electron  emis¬ 
sion  code  was  used  in  the  computations  (Ref.  13). 

The  results  presented  here  are  questionable  in  the  very  low  energy 
limit  (below  several  keV)  because  of  inaccuracies  in  the  QUICKE2  cross 
sections.  T  low-order  Legendre  polynomial  representation  of  angular  dis¬ 
tributions  u  1  in  the  formulation  results  in  anomalies  at  high  energies 
for  low-atomi*.  umber  materials  for  quantities  dependent  on  backward- 
emitted  electrons.  For  example,  the  QUICKE2  code  gives  negative  values 
for  interface  dose  at  the  surface  of  mylar  for  photon  energies  of  600  keV 
and  up.  An  upgraded  QUICKE  code  version  would  solve  the  angular  distribu¬ 
tion  problem  and  extend  cross  section  accuracy  down  to  1  keV  (Ref.  2) .  The 
analytical  developments  of  Strickland  (Ref.  3)  for  electron  energies  below 
1  keV  could  also  be  employed  to  improve  accuracy  in  the  low-energy  range. 

The  dipole  moment  computation  employed  here  for  cable  effects  is  the 
sum  over  energy  of  the  product  of  the  charge  in  each  energy  bin  and  its 
mean  penetration  distance  into  a  dielectric  material.  The  material  used 
here  is  carbon,  although  electron  ranges  are  relatively  insensitive  to 
material  for  low  atomic  numbers  (Ref.  14).  Electrons  are  assumed  to  be 
emitted  normal  to  the  surface  for  this  computation.  A  better  assumption 
would  be  to  use  a  distribution  proportional  to  cos  0  measured  from  the  sur¬ 
face  normal  (Ref.  15).  This  simplification  introduces  factor-of-two  inac¬ 
curacies  in  the  absolute  magnitude  of  the  dipole  moment;  these  inaccuracies 
are  quite  uniform  with  photon  energy  and  should  not  impair  spectrum  sensi¬ 
tivity  comparisons. 

Algorithms  for  calculating  the  excitation  parameters  are  given  in 
Table  5.  Values  for  arbitrary  spectra  can  be  obtained  by  the  expressions 
once  the  terms  for  monoenergetic  spectra  are  known.  These  are  tabulated 
for  representative  configurations  below;  represents  the  energy  fluence 
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Table  5.  SGEMP  Excitation  Term  Algorithms  for  Arbitrary  Spectra 
of  Energy  Fluence  $i  versus  Photon  Energy  i. 
(Lower-case  symbols  represent  the  parameter  value 
computed  for  monoenergetic  spectra  of  unit  fluence.) 


Dose 


D  = 


d. 

1 


Electron  Yield 


Y  = 


V  y- 

1  'x 


£  **  y. 


Average  Velocity  <v>  = 


i  h 


Average  Energy  <  i>  =  — ^ — r — - - 

2-,  i  yi 


Dipole  Moment 


p  =  £  *i  Pj 


of  the  spectrum  in  bin  i.  The  nonlinear  cavity  response  algorithm  is  not 
so  simple,  and  is  presently  in  computer  code  form  (see  Section  5). 

4.1  MONOENERGETIC  SPECTRA 

Excitation  parameters  are  presented  for  monoenergetic  photon  spectra 
in  Appendix  A.,  Results  are  given  in  both  graphic  and  tabular  form.  The 
tables  may  be  used  with  the  algorithms  of  Table  5  to  obtain  excitation 
parameter  values  for  arbitrary  spectra.  A  linear  combination  of  contri¬ 
butions  over  all  energies  in  the  spectrum  using  the  monoenergetic  spectrum 
parameters  as  weighting  factors  gives  the  value  corresponding  to  the  arbi¬ 
trary  spectrum.  An  example  employing  the  tables  is  treated  in  Section  6, 
and  automated  calculations  employing  them  are  described  in  Appendix  D. 
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4.2  BLACKBODY  SPECTRA 


Excitation  parameters  are  given  here  for  blackbody  spectra.  The  spec 
tra  and  parameters  are  displayed  as  functions  of  blackbody  temperature  in 
Appendix  B.  Results  for  arbitrary  combinations  of  blackbodies  may  be 
obtained  by  direct  superposition  of  the  parameter  values  in  the  same  pro¬ 
portion  as  the  fluences  of  the  incident  spectra. 


5.  RESPONSE  PARAMETER  CALCULATIONS 


Satellite  excitation  parameters  have  been  described  above  and  tabulated 
in  the  appendices.  The  results  are  valid  at  all  fluence  levels,  and  are 
dependent  only  upon  attenuation  and  emission  materials  for  a  given  spectrum. 
The  parameters  chosen  to  meter  satellite  response  to  the  excitations  are 
described  here.  These  differ  from  the  excitation  parameters  in  that  they 
require  additional  geometry  specification,  and  they  depend  upon  photon  flux 
rather  than  on  total  fluence.  They  also  may  become  nonlinear  at  higher 
fluence  levels  due  to  space-charge-limiting.  Thus,  the  response  parameters 
are  less  general  than  excitation  parameters,  and  they  are  not  derivable  from 
simple  linear  combinations  of  monoenergetic  results  in  the  space-charge- 
limited  regime. 

A  first-order  measure  of  satellite  response  variations  due  to  varying 
photon  spectrum  may  be  achieved  with  a  self-consistent  solution  of  electron 
emission  between  grounded  parallel  plates.  This  solution  resembles  inter¬ 
ior  cavity  response  while  remaining  general  and  relatively  simple  to  compute. 

2 

In  addition,  results  scale  with  the  parameter  Fd  ,  where  F  is  the  photon  flux 
and  d  is  the  diode  spacing.  Results  are  generalized  in  that  they  apply  to  a 
continuous  range  of  spacings  rather  than  to  one  specific  case  (see  Ref.  1 J . 

Response  parameters  employed  here  are  the  electric  field  at  the  front 
plate,  the  peak  cavity  potential,  the  transmitted  current,  and  the  electron 
charge  dipole  moment.  The  dipole  moment  presented  here  is  the  integral  over 
the  entire  charge  cloud  of  the  product  of  the  charge  density  times  its  dis¬ 
tance  from  the  emission  face.  Electron  energy  spectra  are  obtained  from  the 
QUICKE2  code.  Angular  variation  of  the  electron  number  intensity  is  propor¬ 
tional  to  cos  0  where  0  is  measured  from  the  surface  normal.  Charge  is 
emitted  from  the  forward  plate.  The  treatment  is  equivalent  to  that  in 
Reference  16,  extended  to  arbitrary  emission  energy  spectra  between  paral- 
el  plates. 

Results  for  blackbody  and  blackbody  combination  spectra  incident  on 
aluminum  plates  are  shown  in  Appendix  C.  The  combination  spectra  plotted 
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were  chosen  from  a  larger  set  of  combinations  in  the  range  2  to  35  keV; 
the  ones  chosen  bracket  the  responses  at  virtually  all  levels. 

Magnetic  field  estimates  can  be  immediately  made  from  the  results 
using  the  simple  expression: 


where 

=  magnetic  field  (amp/m)  , 

J  =  transmitted  current  density  (amp/m^)  , 
R  =  diode  radius  (m)  . 


6.  EXAMPLE  USES  OF  THE  PARAMETERS 


Representative  satellite  excitation  parameters  for  SGEMP  have  been 
determined  and  evaluated  in  previous  sections  of  this  report.  Example  uses 
of  the  parameters  are  discussed  here.  Excitation  and  response  parameters 
for  several  bremsstrahlung  photon  spectra  are  computed  and  compared  with 
values  for  a  nominal  spectrum.  First,  a  sample  illustration  of  the  numer¬ 
ics  is  given. 

6.1  SAMPLE  COMPUTATION 

A  sample  computation  illustrating  the  use  of  the  monoenergetic  excita¬ 
tion  parameters  is  shown  in  Figure  1.  The  electron  yield  from  fiberglass 
behind  a  solar  array  is  obtained  for  the  spectrum  shown  by  dividing  it  into 
energy  bins  of  width  20  keV  (the  large  bin  width  is  used  for  illustrative 
purposes  only,  to  reduce  the  numerical  computation  effort) .  The  fluence  in 
each  energy  bin  and  the  electron  yield  corresponding  to  that  energy  are 
obtained  (indicated  by  dots  on  the  curves) .  The  product  of  the  two  gives 
the  electron  yield  due  to  photons  in  that  energy  bin.  The  sum  over  all 
bins  gives  the  electron  yield  for  the  entire  spectrum.  These  steps  are 
illustrated  mathematically  by  the  relation: 

N 

Y  =  y'  $  Ac.  y. 

Z-f  e,  1  'l 

i=l 

*  20  keV  x  (2.7xl0-2  x  2.0xl0"8  +  1.6xl0-2  x  1.2xl0“8  +  4.5xl0"3 

x  6.8xl0-9  +  8.7xl0"4  x  4.5xl0-9) 

s  1.5  x  10"8  C/cal  (QUICKE  ...  1.4  x  10-8  C/cal)  . 

In  this  particular  case,  the  results  obtained  with  the  large-energy  bins 
correspond  very  well  to  results  obtained  with  QUICKE2  using  much  finer  bins. 
This  approach  applies  to  every  parameter  type  using  the  corresponding  algo¬ 
rithms  of  Table  5.  Notice  that  we  have  assumed  unit  fluence  to  obtain  the 
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Figure  1.  Sample  calculation  of  electron  yield  due  to  an  8-keV 
blackbody  incident  on  a  solar  array.  Monoenergetic 
photon  results  are  employed. 


units  C/cal  for  purposes  of  comparing  spectra  (see  Figure  1).  This  nor¬ 
malization  is  employed  for  all  excitation  parameter  results  in  this  report. 

6.2  BREMSSTRAHLUNG  SPECTRUM  EVALUATIONS 

A  subset  of  the  parameters  defined  in  previous  sections  has  been  eval¬ 
uated  for  bremsstrahlung  photon  spectra  in  the  range  from  135  to  400  keV. 

The  spectra  are  shown  in  Figure  2.  The  parameters  used  are.  indicated  in 
Figure  3,  below  the  spectrum  profiles  presented  as  ratios  of  the  corres¬ 
ponding  values  for  a  10-keV  blackbody.  The  profiles  diverge  from  unity  as 
the  bremsstrahlung  energy  increases.  Examination  of  the  10-keV  blackbody 
spectrum  (see  Appendix  B)  and  the  bremsstrahlungs  corroborates  this  behav¬ 
ior  .  The  lower-energy  bremsstrahlung  spectra  are  more  similar  in  energy 
content  to  the  10-keV  blackbody  than  the  higher-energy  bremsstrahlungs. 

The  cavity  response  to  the  incident  bremsstrahlungs  as  a  function  of 
incident  flux  is  shown  in  Figures  4  and  5.  The  spectra  penetrate  3  mils  of 
aluminum,  emitting  electrons  to  the  interior  from  the  first  plate.  Some 
space-charge-limiting  is  evident  from  the  curve  nonlinearities.  For  a  diode 

spacing  of  0.5  m  and  a  pulse  width  of  30  nsec,  space-charge-limiting  begins 

-4  2 

at  about  5  x  10  cal/ cm  . 


RATIO  BREMSSTRAHLUNG/lO-keV  BLACKBODY 


I 

i 


$ 


RE-01451 


RESPONSE  PARAMETERS 


1  YIELD,  3-mil  AS. 

2  YIELD,  10-mil  MYLAR 

3  YIELD,  SOLAR  ARRAY 

4  <v>  BEHIND  3-mil  Ail 

5  YIELD,  FORE/BACK  BEHIND  SOLAR  ARRAY 

6  YIELD,  HI/LO-Z  BEHIND  SOLAR  ARRAY 


7  YIELD,  20-mil  AS,  BOX 

8  DOSE  INSIDE  DEVICE  HOUSING 

9  DOSE  RATIO  Au/Si  INSIDE  HOUSING 

10  DOSE  IN  SHIELDED  BOX 

11  CABLE  YIELD  (3  mils) 

12  CABLE  DIPOLE  MOMENT  (3  mils) 


Figure  3.  Ratios  of  bremsstrahlung  spectrum  excitation  parameters  to  corres¬ 
ponding  values  for  a  10-keV  blackbody.  The  bremsstrahlung  voltage 
V  is  indicated  for  each  figure. 
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Figure  4.  Diode  potential  and  scaled  current  due  to  incident  bremsstrahlung 
spectra  of  voltage  V  (keV) .  F  is  the  photon  flux,  J  is  the  elec¬ 
tron  current  density,  and  d  is  the  plate  spacing. 
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Figure  5.  Dipole  moment  and  scaled  electric  field  at  the  front  face  of  a 
diode  due  to  incident  bremsstrahlung  spectra  of  voltage  V  (keV) 
F  is  the  photon  flux,  E  is  the  electric  field,  and  d  is  the 
plate  spacing. 
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Figure  A-l.  Interface  dose  at  the  surface  of  dielectric  materials 
due  to  monoenergetic  photons  incident  at  1  cal/cm 
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Figure  A-2.  Reverse-emitted  electron  yield  at  the  surface  of  quartz 
due  to  monoenergetic  photons  incident  at  1  cal/cm2 
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Figure  A-3.  Interface  dose  behind  thermal  blanket  due  to  monoener 
getic  photons  incident  at  1  cal/cm2 
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Figure  A-4.  Electron  yield  behind  thermal  blanket  due  to  monoener- 
getic  photons  incident  at  1  cal/cm2 
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Figure  A-5.  Forward-emitted  electron  yield  from  aluminum  after  pene¬ 
tration  of  3  and  20  mils  aluminum  by  moroenergetic  photons 
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Figure  A-6.  Average  velocity  and  energy  of  electrons  emitted  from 
aluminum  due  to  monoenergetic  photons  incident 
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Figure  A-7.  Electron  yields  from  fiberglass  behind  solar  array  due 
to  monoenergetic  photons  incident  at  1  cal/cm2 
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Figure  A-8.  Reverse- emitted  electron  yields  from  gold  and  fiberglass 
behind  solar  array  due  to  monoenergetic  photons  incident 
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Figure  A-10.  Electron  yield  inside  cables  of  shield  thicknesses  3 
and  20  mils  due  to  monoenergetic  photons  incident  at 
1  cal/cm2 
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Figure  A-ll.  Dipole  moment  inside  cable  dielectric  due  to  monoenergetic 
photons  incident  at  1  cal/cm2.  Copper  shield  thicknesses 
are  indicated  beside  the  curves. 
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Figure  A-12.  Bulk  dose  inside  box  due  to  monoenergetic  photons  inci- 
i  I  dent  at  1  cal/cm2 
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Figure  A-13.;  Bulk  dose  ratio  in  gold  compared  to  silicon  for  mono 
energetic  photons 
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Figure  A-14.  Electron  yields  forward- emitted  from  materials  inside  box 
due  to  monoenergetic  photons  incident  at  1  cal/cm2 


£ 


43 


RE-01442 


™0N  ELECTRON  YIELD  (C/cal) 


ceV) 

FORWARD 

REVERSE 

2 

_ 

_ 

5 

_ 

_ 

10 

1.797xl0"8 

1.359x10'° 

15 

9.154x10"° 

6.381x10"° 

20 

9.410xl0"8 

6.159x10"° 

30 

5.408x10"° 

3.172x10"° 

40 

3.060x10"° 

1.655x10"° 

50 

2.020x10'° 

9.382xl0"9 

60 

1.360x10"° 

6 . 1 20x1 O-9 

70 

9.930xl0'9 

4 . 09 1 xl 0"9 

80 

7.708xl0"9 

2.905xl0-9 

PHOTON 

ENERGY 

(keV) 


ELECTRON  YIELD  (C/cal) 


FORWARD 


REVERSE 


90 

100 

125 

150 

200 

300 

400 

500 

600 

800 

1000 


6.372x10 

5.542x10 

4.835x10 

4.835x10 

5.967x10 

9.062x10 

1.202x10 

1.451x10 

1.653x10 

1.959x10 

2.155x10 


-9 

-9 

-9 

-9 

-9 

-9 

-8 

-8 

-8 

-8 

-8 


2.178x10 

1.724x10 

1.204x10' 

1.015x10' 

1.065x10 

1.495x10' 

1.922x10 

2.271x10 

2.534x10' 

2.864x10 

3.126x10 


-9 

-9 


-9 


-9 

-9 


-9 

-9 


Figure  A-15.  Electron  emission  from  aluminum  behind  box  wall  due  to 
monoenergetic  photons  incident  at  1  cal/cm2 
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Figure  A-16.  Electron  emission  ratios  due  to 

monoenergetic  photons. 

The  high-to-low-Z  ratio  is  for  forward  emission  from 
gold  compared  to  carbon.  The  forward -to-reverse  ratio 
is  for  aluminum. 
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Figure  A-17.  Bulk  dose  inside 
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photons  incident  at  1  cal/cm2.  The  ratio  of  the  two 
curves  is  shown  in  Figure  A-13. 


46 


Z  M 


RE-01224 


PHOTON  ENERGY  (keVl 


PHOTON 

ENERGY 

(keV) 

BULK  DOSE  ? 
[rad(Si  )/cal/cnr] 

PHOTON 

ENERGY 

(keV) 

BULK  DOSE  ? 
[rad(Si )/cal/cnr] 

2 

_ 

90 

1.202x10° 

5 

- 

100 

1.980x10s 

10 

- 

125 

3.985xlOS 

15 

- 

150 

5.620xl0S 

20 

- 

200 

7, 904xl03 

30 

1.332x10° 

300 

1 .01 lxl 04 

40 

5.340X102 

400 

1.1 00x1 0^ 

50 

3.771xl03 

500 

1.1 58x1 04 

ft 

60 

6.607xl03 

600 

1. 205x1 0H 

70 

1 . 705x1 03 

800 

1.209xl04 

80 

6.973xl02 

1000 

1.1 69x1 04 

Figure  A-18.  Dose  inside  shielded  box  due  to  monoenergetic  photons 
incident  at  1  cal/cm2 
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APPENDIX  B 


EXCITATION  PARAMETER  RESULTS 
FOR  BLACKBODY  SPECTRA 
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Figure  B-2.  Interface  dose  at  the  surface  of  dielectric  materials 
due  to  blackbody  photon  spectra  incident  at  1  cal/cm2 
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Figure  B-3.  Electron  yield  at  the  surface  of  quartz  due  to  black 
body  photon  spectra  incident  at  1  cal/cm2 
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figure  B-4.  Interface  dose  behind  thermal  blanket  due  to  blackbody 
photon  spectra  incident  at  1  cal/cm2 
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Figure  B-5.  Electron  yield  behind  thermal  blanket  due  to  black 
body  photon  spectra  incident  at  1  cal/cm2 
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Figure  B-6.  Electron  yield  forward- emitted  from  aluminum  after 
penetration  of  3  and  20  mils  aluminum  by  blackbody 
photon  spectra 
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Figure  B-7.  Average  velocity  and  energy  of  electrons  emitted  from 

3  mils  aluminum  due  to  blackbody  photon  spectra  incident 
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Figure  B-8.  Electron  yield  behind  solar  array  due  to  blackbody 
photon  spectra  incident  at  1  cal/cm2 
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Figure  B-9.  Electron  emission  ratios  behind  solar  array  due  to 
blackbody  photon  spectra  incident 
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Figure  B-10. 


Electron  yield  inside  cables  of  shield  thicknesses 
3  and  20  mils  due  to  blackbody  photon  spectra  inci 
dent  at  1  cal/cm2 
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Figure  B-ll.  Dipole  iroment  inside  cable  dielectric  due  to  blackbody 
photon  spectra  incident  at  1  cal/cm2.  Copper  shield 
thicknesses  are  indicated  beside  the  curves. 
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Figure  B-12.  Bulk  dose  and  dose  ratio  inside  box  due  to  blackbody 
photon  spectra  incident  at  1  cal/cm2 
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Figure  B-13.  Forward -emit ted  electron  yield  inside  box  duetto 
blackbody  photon  spectra  incident  at  1  cal/cm 


Figure  B-14.  Electron  emisson  ratios  inside  box  due  to  blackbody 
photon  spectra.:  The  high-to-low-Z  emission  ratio  is 
for  forward  emission  from  gold  compared  to  carbon. 
The  forward-to- reverse  ratio  is  for  aluminum. 
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Figure  B-15.  Bulk  dose  and  dose  ratio  inside  device  housing  due 
to  blackbody  photon  spectra  incident  at  1  cal/cm2 


64 


<->  10 


1  cal/cnr 


ce 

c 

> 

^  o 
<  ^ 


</) 

UJ 

r— 

r— 

r— 

CO 

•n 

•r— 

■  r— 

o 

£ 

£ 

s 

Q 

cm 

o 

o 

•r— 

r— 

CM 

CO 

RE-01213 


10  20 
BLACKBODY  TEMPERATURE  (keV) 


BLACKBODY 

TEMPERATURE 

(keV) 

DOSE 

[rad(Si)] 

5 

1.34xl02 

8 

8.93xl02 

10 

1 .41xl03 

15 

2.1 9x1 03 

20 

2.69xl03 

Figure  B-16.  Dose  inside  shielded  box  due  to  blackbody  photon 
spectra  incident  at  1  cal/cm2 
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RESPONSE  PARAMETER  RESULTS 
FOR  BLACKBODY  SPECTRA 
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Figure  C-l.  Diode  potential  and  scaled  current  due  to  incident  blackbody 
spectra  of  temperature  T  (keV) .  F  is  the  photon  flux,  J  is 
the  electron  current  density,  and  d  is  the  plate  spacing. 
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Figure  C-2.  Dipole  moment  and  scaled  electric  field  at  the  front  face  of  a 
diode  due  to  incident  blackbody  spectra  of  temperature  T  (keV) 
F  is  the  photon  flux,  E  is  the  electric  field,  and  d  is  the 
plate  spacing. 
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Figure  C-3.  Diode  potential  and  scaled  current  due  to  blackbody  combination 
spectra.  Equal  fluxes  of  temperature  T  (indicated  in  keV)  are 
incident.  F  is  the  photon  flux,  J  is  the  electron  current,  and 
d  is  the  plate  spacing. 
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Figure  C- 
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4.  Dipole  moment  and  scaled  electric  field  at  the  front  face  of  a 
diode  due  to  blackbody  combination  spectra.  Equal  fluxes  of 
temperature  T  (indicated  in  keV)  are  incident.  F  is  the  photon 
flux,  E  is  the  electric  field,  and  d  is  the  plate  spacing. 
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The  excitation  parameters  discussed  in  Section  4  of  this  report  are 
conveniently  obtainable  for  arbitrary  photon  spectra  by  means  of  program¬ 
mable  calculator  codes  described  here.  The  programs  are  designed  for  use 
on  the  Texas  Instruments  TI-59  calculator  and  are  available  on  magnetic 
cards.  Arbitrary  spectra  can  be  keyed  in  or  read  in  from  magnetic  cards, 
and  the  calculator  will  automatically  compute  the  excitation  parameters. 
Brief  descriptions  of  the  codes  and  specific  user  instructions  follow. 

D.l  CODE  OVERVIEW 

A  brief  overview  of  the  codes  is  helpful  in  understanding  their  oper¬ 
ation.  The  calculator  memory  is  divided  into  four  banks  which  can  be  pro¬ 
grammed  with  up  to  240  steps  each.  Memory  can  be  partitioned  to  interchange 
code  steps  and  data  storage,  with  up  to  30  numbers  stored  in  each  bank.  One 
insertion  of  a  magnetic  card  sets  all  locations  in  a  bank  of  memory.  There¬ 
fore,  the  memory  has  been  programmed  according  to  the  divisions  shown  in 
Table  D-l.  The  memory  allocation  shown  limits  the  program  to  only  240 
steps,  but  it  has  the  considerable  advantage  of  permitting  spectra  to  be 
stored  on  one  set  of  cards  and  excitation  parameter  data  on  another  set. 

The  user  may  read  in  a  spectrum  to  bank  2  and  compute  parameters  for  differ¬ 
ent  satellite  configurations  by  simply  reading  in  different  data  to  bank  3 
and  repeating  the  calculation. 


Table  D-l.  TI-59  Memory  Allocation  for  Computation  of 
Spectrum  Evaluation  Parameters 


Bank 

1 

Program 

Up  to 

240  steps 

Bank 

2 

Spectrum 

Up  to 

14  point  pairs 

Bank 

3 

Excitation 
parameter  curves 

Up  to 

15  point  pairs 

Bank 

4 

Working  storage 

D.2  USER  INSTRUCTIONS 

Step-by-step  instructions  for  using  the  programs  are  given  below.  The 
user  is  assumed  to  have  an  elementary  knowledge  of  the  calculator. 
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Partition  calculator  memory  to  240  maximum  program  steps  with  90  data 
locations : 

•  Press  9  Op  17 

•  Press  CLR 

•  Insert  program  card  (called  ARBSPC) 

Calculator  will  display  1,  indicating  bank  1  has  been  programmed. 


D.2.2  Spectrum  Input 

•  Punch  spectrum  into  memory  bank  2: 

Energies  (keV)  locations  60-73 

Differential  fluence  (energy/ energy/ area)  locations  75-88 

(To  save  thi«  spectrum  on  a  card,  press  CLR,  then  WRITE,  and 
insert  a  blank  card.) 

Up  to  14  energies  and  14  differential  fluences  may  be  input.  A  spec¬ 
trum  may  be  read  in  from  a  card  by  pressing  CLR  and  inserting  the  card.  The 
calculator  will  display  2.,  indicating  memory  bank  2  has  been  programmed. 

See  optional  input  routine,  TLBSET,  for  rapid  spectrum  input.  Spectra  are 
normalized  automatically  by  the  code.  Photon  energies  must  be  less  than 
1000  keV.  Spectra  must  have  at  least  three  energy  bins. 

D.2.3  Excitation  Parameter  Data  Input 

Read  selected  excitation  parameter  weighting  factor  (see  Appendix  A) 
into  locations  30-59  (bank  3).:. 

•  Press  CLR 

•  Insert  card 

Calculator  displays  3.,  indicating  bank  3  has  been  programmed. 

D.2.4  Execution 

The  calculator  is  now  ready  to  compute  the  selected  parameter  for  the 
spectrum  chosen. 

•  Press  RST 

•  Press  R/S 


Execution  commences  and  numbers  flash  on  the  display  for  each  energy  bin  in 
the  order  below. 

ith  bin  energy  (keV) 

$Qe  ith  differential  fluence  (energy/ energy/area) 

53  $q  Ae^  cumulative  sum  of  fluences  in  each  bin  (energy/ area) 
i  ei 

When  the  display  becomes  fixed,  the  value  is  total  fluence: 


Nbins 

Ac.  (energy/area) 
i=l  ei  1 

where  Nbins  is  the  total  number  of  energy  bins.  This  is  the  total  fluence 
incident  on  the  exterior  of  the  excitation  parameter  configuration  chosen. 
Press  R/S 

Calculator  displays  Excitation  Parameter  P  (per  unit  fluence) . 

Press  R/S 

Calculator  displays  ratio  Pprevious/pcurrent>  the  ratio  of  the  previously 
calculated  excitation  parameter  to  the  current  one.  Thus,  if  a  dose  ratio 
Au/Si  is  being  computed,  for  example,  the  dose  for  gold  is  computed  first, 
then  a  new  card  read  in  for  silicon  and  the  process  repeated.  The  third 
static  display  yields  the  ratio  on  the  second  pass. 

The  execution  steps  are  summarized  in  Table  D-2. 


D.2.5 


itional  Spectrum  Input  Routine 


A  routine  is  available  to  expedite  table  entries  and  checkouts.  In 
our  case,  we  want  to  input  a  string  of  energies  into  memory  beginning  at 
location  60,  and  the  corresponding  differential  fluences  beginning  at  loca¬ 
tion  75. 

To  perform  this  operation,  press  9  Op  17  to  partition  memory,  then 
read  in  the  TLBSET  code  card.  Press  RST,  then  R/S.  The  calculator  now 
displays  1 . ,  indicating  energy  bin  1 .  Energies  may  now  be  input  by  keying 
them  in  and  pressing  R/S.  The  index  of  the  next  bin  then  appears.  To 
indicate  the  last  bin,  enter  a  zero  and  the  calculator  will  automatically 
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Table  D-2.  Summary  of  Program  Setup  and  Execution,  ARBSPC 


S  *0e.  Aei 

1  1 

Calculator  stops  and  displays: 

Nbins 

X!  Sp,  Ae.  (total  fluence  incident) 
i=l  1  1 

2 

Press  R/S.  Calculator  displays  P  (excitation  parameter/cal/cm  ) 

Press  R/S.  Calculator  displays  ratio  P  _  .  J P  „ 

— —  r  previous  current 

To  continue. 

Press  CLR 

Read  in  new  spectrum  or  data  table 
Press  RST 
Press  R/S 


skip  to  location  75  and  display  1.  for  bin  1.  At  this  point,  begin  entering 
differential  fluences. 

The  values  may  be  displayed  in  succession  by  pressing  E  and  then  R/S, 
R/S,  etc.  Displays  beginning  at  any  location  can  be  commenced  by  entering 
the  desired  location  and  pressing  A'.  The  calculator  will  automatically 
begin  displaying  at  location  60  if  no  entry  is  made  into  A'. 

A  summary  of  TLBSET  operations  is  given  in  Table  D-3.  A  summary  of 
calculator  outputs  is  given  in  Table  D-4. 


Table  D-3.  Summary  of  Input/Output  Routine  TLBSET  Instructions 


k' 

B  ' 

C' 

D' 

E' 

Starting 
location  for 
display 
(default  = 
value  in  B) 

Starting 
location  to 
multiply,  R/S, 
multiply  by, 
R/S 

B 

C 

D 

E 

Starting 
location  for 
storage 
(default=60) 

Display 
values 
beginning 
at  A' 

•  Press  9  Op  17 

•  Read  in  program  TLBSET 

•  Press  RST,  R/S 

•  Enter  beginning  storage  location  in  (default  value  =  60) 

•  Enter  energy,  press  R/S,  energy,  R/S,  etc. 

•  Terminate  string  with  0^,  R/S 

•  Enter  differential  fluence,  press  R/S,  etc. 

•  To  display  values,  press  E,  R/S,  R/S,  .  .  . 

To  display  values  beginning  at  a  location  other  than  B,  enter 
the  location  in  k'  after  RST,  then  press  E,  R/S,  R/S,  .  .  . 
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Table  D-4.  Spectrum  Evaluation  Program  Symbols  and  Units 


Energy 

t 

keV 

Photon  fluence 

$ 

cal/cm^ 

Photon  differential  fluence 

$  £ 

2 

cal/cm  /keV 

Dose 

D 

rad (material) 

Electron  yield 

Y 

C/cal 

Velocity 

V 

m/sec 

Electron  dipole  moment 

P 

C-m/cal 
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